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ABSTRACT 

Integral field spectroscopy of the LV2 proplyd is presented taken with the 
VLT/FLAMES Argus array at an angular resolution of 0.31x0.31 arcsec^ and ve- 
locity resolutions down to 2 kms~^ per pixel. Following subtraction of the local M42 
emission, the spectrum of LV 2 is isolated from the surrounding nebula. We measured 
the heliocentric velocities and widths of a number of lines detected in the intrinsic 
spectrum of the proplyd, as well as in the adjacent Orion nebula falling within a 6.6 
X 4.2 arcsec^ field of view. It is found that far-UV to optical collisional lines with crit- 
ical densities, iVcr, ranging from lO'^ to 10^ cm~^ suffer collisional de-excitation near 
the rest velocity of the proplyd correlating tightly with their critical densities. Lines 
of low A^cr are suppressed the most. The bipolar jet arising from LV2 is spectrally 
and spatially well-detected in several emission lines. We compute the [O ill] electron 
temperature profile across LV 2 in velocity space and measure steep temperature vari- 
ations associated with the red-shifted lobe of the jet, possibly being due to a shock 
discontinuity. From the velocity-resolved analysis the ionized gas near the rest frame 
of LV2 has To = 9200 ± 800 K and N^ - lO^cm-^, while the red-shifted jet lobe 
has Te ss 9000 - lO^K and Ac - 10^ - lO'^cm-^. The jet flow is highly ionized but 
contains dense semi-neutral clumps emitting neutral oxygen lines. The abundances of 
N+, 0^+, Ne^+, Fe^+, S+, and S^+are measured for the strong red-shifted jet lobe. 
Iron in the core of LV2 is depleted by 2.54 dex with respect to solar as a result of 
sedimentation on dust, whereas the efficient destruction of dust grains in the fast mi- 
crojet raises its Fe abundance to at least 30 per cent solar. Sulphur does not show 
evidence of significant depletion on dust, but its abundance both in the core and the 
jet is only about half solar. 

Key words: ISM - abundances; HII regions; ISM: individual objects - (LV2, 167-317, 
Orion Nebula); stars: pre-main-sequence; protostars; planets and satellites: protoplan- 
etary disks 



1 INTRODUCTION 

Protoplanetary disks (proplyds) provided the first evidence 
of gaseous dusty disks around young stars in the early 1990s. 
The archetypal proplyds found in the Orion nebula (M42; 
O'Dell et al. 1993) are dense semi-ionized objects photoevap- 
orated to various degrees depending on their distance from 
the ionizing Trapezium cluster. They are landmark objects 



"^ Based on observations made with ESO telescopes at 
Paranal Observatory under programme 078.C-0247(A). 
f E-mail: ytsamis@eso.org 
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in the study of how circumstellar disks and eventually plan- 
etary systems form (e.g. O'Dell 2001). They are a unique 
environment for the study of disk formation and evolution 
in areas dominated by massive OB-type stars. Massive stel- 
lar clusters and their associated H ll regions, such as in the 
Orion OBI association, are thought to represent the clos- 
est analogs to the birth environment of our Solar System 
(Adams 2010). The proplyds of Orion should therefore be 
ideal laboratories for the study of our planetary system's 
origins. At optical wavelengths proplyds present a photoion- 
ized skin facing the ionizing cluster, giving way to a dusty 
envelope often shaped into comet-shaped outflows. The el- 
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Table 1. Journal of VLT Argus observations.' 



Date 


A-range 


Grating 


A/5A 


Exp. time 


(UT) 


(A) 






(sec) 


2006/10/08 


6383-6626 


H14B 


46 000 


3x139 


2006/12/30 


3700-3867 


HRl 


36 000 


3x188 


2006/10/10 


4033-4201 


HR3 


39 000 


3x185 


2006/12/30 


4033-4201 


HR3 


39 000 


3x185 


2007/01/09 


4188-4392 


HR4 


32 500 


3x296 


2007/01/09 


4538-4759 


HR6 


32 500 


3x225 


2007/01/09 


4917-5163 


HR8 


32 000 


3x130 



" The Argus array was centered at (RA, Dec)jD2000 = 
(05''35"'16.857'', -05°23'15.03") at a position angle of 



emental content and chemistry of proplyds are very poorly 
known, however, (i) the composition of planet-forming cir- 
cumstellar envelopes is of great topical interest given the es- 
tablished positive correlation between host star metallicity 
and the incidence of giant planetary companions (e.g. Gon- 
zalez 1997; Neves et al. 2009); (ii) higher metallicity in the 
protoplanetary disk favours the formation of giant planets in 
the 'core accretion' scenario (e.g. Boss 2010). While the pre- 
viously discrepant abundances for the interstellar medium 
and main-sequence stars in Orion are beginning to converge 
(Simon-Diaz & Stasinska 2011), proplyds are the last major 
component of Orion for which a chemical abundance scale 
is lacking. 

We have been taking steps to rectify this situation with 
a programme targeting a sample of bright proplyds in M42. 
In Tsamis et al. (2011; Paper I hereafter) a chemical abun- 
dance study of the protoplanetary disk LV 2 (Laques & Vi- 
dal 1979) and its Orion nebula host vicinity was presented, 
based on the analysis of VLT optical integral field spec- 
troscopy and HST FOS single aperture ultraviolet to far-red 
spectroscopy. From an emission-line analysis the abundances 
of several elements were measured for the proplyd and the 
local M42 nebular field. LV 2 was found to be slightly over- 
abundant in carbon, oxygen and neon compared to the Orion 
nebula gas-phase composition and to the Sun. The carbon, 
oxygen and neon abundances in LV2 were measured to be 
«0.2-0.3 dex higher than those in B-type stars of the Ori 
OBI association studied by Simon-Diaz (2010) and Simon- 
Diaz & Stasinska (2011). 

That result constitutes a direct measure of the metallic- 
ity of gas photo-evaporated from circumstellar disk material 
where planet formation may be underway. 

In this paper, we present higher velocity dispersion inte- 
gral field spectra of LV 2 which enable a physical analysis of 
its fast microjet to be undertaken. The paper is organized as 
follows. The dataset specifications and reduction method are 
outlined in Section 2. The analysis and results are presented 
in Section 3, followed by our conclusions in Section 4. 



2 OBSERVATIONS AND REDUCTIONS 

Integral field spectroscopy of LV 2 was performed on the 8.2- 
m VLT/UT2 Kueyen during 2006 October and December, 
and 2007 January with the FLAMES Giraffe Argus array. A 
field of view of 6.6 x 4.2 arcsec^ was used yielding 297 posi- 



tional spectra in the optical range from six high dispersion 
grating settings (Table 1). The size of the angular resolution 
element was 0.31 x 0.31 arcsec^, corresponding to a spatial 
scale of 123 x 123 AU^ at the distance to M42 (412 pc; 
Reid et al 2009) . The typical seeing during the observations 
was <0.8 arcsec FWHM. The data were cosmic ray cleaned, 
flat fielded, wavelength calibrated, and extracted with the 
girBLDRS pipeline also employed in Paper I (see Blecha & 
Simond 2004 for details) . The resulting data cubes have not 
been corrected for the effects of differential atmospheric re- 
fraction (DAR; as in Paper I), as for these high dispersion 
spectra the wavelength coverage per grating is small and the 
effects of DAR are negligible under the low (<1.5) observed 
air- masses. 

The flux calibration was done within iraf using expo- 
sures of various spectrophotometric standards for the vari- 
ous high resolution (HR) gratings: CD-329927 (HRl), LTT 
1020 (HR3, HR8), LTT 7987 (HR14B) - Hamuy et al. 
(1992), Hamuy et al. (1994), Feige 67 (HR4) - Oke (1990). 
These standards are accurate to about one per cent for rela- 
tive flux calibration and a few per cent in absolute flux units 
- under the non-photometric observing conditions, however, 
the absolute flux calibration cannot be more accurate than 
about flve per cent. The transfer of relative flux between 
different wavelength ranges adds most to the uncertainties, 
therefore five per cent accuracy may be a conservative esti- 
mate. 

The FWHM velocity resolution of the spectra derived 
from measurements of comparison Th-Ar arc lines was 7.68 
± 0.10 kms"^ at 656.3 nm (HR14B grating), 10.1 ± 0.1 
kms"^ at 372.6 nm (HRl), 8.60 ± 0.20 kms"^ at 409.0 nm 
(HR3), 9.60 ± 0.10 kms"^ at 434.8 nm (HR4), 9.80 ± 0.10 



kms'^ at 465.8 nm (HR6), and 10.6 ± 0.1 kms" 
nm (HR8). 



at 496.5 



3 RESULTS 

3.1 Extracted spectra and LV2 versus M42 
surface brightness 

The flux calibrated high dispersion Ha spectrum of LV 2 is 
shown in Fig. 1 in the heliocentric velocity frame. The Orion 
nebula contribution to the line flux has not been subtracted; 
the intrinsic line flux from LV 2 is about 65 per cent of the 
total. This image was formed by integrating the signal of 
the {X, Y, A) HR14 data cube along rows X= 2 - 4 of 
the integral fleld array which contain most of the emission 
from the proplyd, and then extracting the ( — 220, -1-260) km 
s~^ velocity range along columns Y — 1 - 22. The proplyd 
peak emission is concentrated along those columns where 
additional emission is detected at around —100 and -1-150 
kms~^, respectively originating in the approaching and re- 
ceding lobes of LV2's jet; these are well-separated from the 
line systemic velocity. 

Extracted ID spectra of LV 2 and of the local Orion 
nebula are shown in Figs. 2 and 3. These were formed by 
co-adding the spectra from nine spaxels of the reconstructed 
data cube containing the peak of LV 2's emission, and about 
70 spaxels where proplyd emission was not detected (based 
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is logarithmic with a miniirLum of 5.4xl0~ ' and a maximum of 1.1x10" "^ erg s~^ cm~^. The length of the vertical axis is 6.6 arcsec 
(22 spaxels). The vertical spaxel size has been magnified by a factor of 2.7 for display purposes. The Orion nebula emission has not been 
subtracted. See the text for details. 



on a Ha image of the field extracted from the cube) [j The 
summed spectra were scaled by the number of spaxels in 
each and the Orion spectrum was subtracted from the for- 
mer spectrum to obtain pure LV2 spectra for the various 
grating settings. The local Orion nebula extraction was at- 
tempted with two different definitions of the background, in 
terms of number of co-added spaxels and their location, and 
the resulting LV 2 spectral line profiles and continuum level 
were found to be very similar. This simple way of perform- 
ing the background subtraction and obtaining the intrinsic 
proplyd emission does not take into account the fact that 
LV 2 may be more dusty than the surrounding nebula. 

The issue of internal extinction in a proplyd has been 
discussed in detail by Henney & O'Dell (1999), who showed 
that neglecting the attenuation of flux due to dust withm 
a proplyd can lead to over-subtraction of the background, 
in the limiting case that the principle emitting layer of the 
Orion nebula lies behind a given proplyd. This can result 
in artificially distorted line profiles to be obtained, for ex- 
ample, with double peaks and negative troughs. The loga- 
rithmic reddening at H/3 within LV 2 was estimated in Pa- 
per I to be ~0.5 dex higher than for the local Orion neb- 
ula, which translates to an effective extinction, Av, about 
1 mag larger. In this data set we do not deduce, however, a 
discernible wavelength-dependent effect on the resulting line 
profiles that could be caused by an apparent reddening ex- 
cess. For example, the high order Hi lines (transitions from 
upper levels n > 10) in the violet part of LV2's spectrum 
(<3800 A) are all similarly shaped and well-fitted by Gaus- 
sians (of FWHM ^ 0.6 A) with no evidence of anomalous 



^ The typical background region was defined as {X, Y) = (4- 
11, 18-20) + (7-12, 4-5) + (11-12, 6-17) and slight variations 
thereof were used to test its effect on the background-subtracted 
LV2 spectra. 



profiles. At the same time, the [O ll] violet doublet lines 
at 3726, 3729 A appear heavily suppressed in the subtracted 
spectrum. Very weak lines at longer wavelengths such as the 
recombination lines C ll A4267 and O ll A4649 (from recom- 
bining C^^ and O^""", respectively), are well-recovered in the 
intrinsic proplyd spectrum (see also Fig 4). The [N ll] A6584 
line is also heavily suppressed, exhibiting a central trough, 
a characteristic not shared by the [S ll] A4069 line. In fact, 
the only lines that show any evidence of altered profile shape 
following the background subtraction are the forbidden lines 
[O ll] AA3726, 3729, [N ll] 6548, 6584, [Ar iv] 4711, 4740, and 
perhaps only very slightly, [O ill] A4959 (Fig 4). 

From this we draw the conclusion that the dominant 
effect on those line profiles that appear distorted in this 
data set is not due to the poor definition or over-subtraction 
of the Orion nebula's local emission, but that colhsionally- 
excited lines (CLs) are quenched near the systemic velocity 
of LV2, due to the high plasma density there. The degree 
of suppression correlates very well with the critical densities 
of the upper levels of the lines in question. In Fig. 5 (left) 
the observed linear relationship is shown between the line 
critical density versus the fraction. Fa, of specific line flux 
emitted from LV2 over the total (LV2 + M42), normalized 
per spaxel. A linear fit was performed through the VLT data, 
excluding for consistency the coloured data points which 
are from the far-UV HST FOS measurements presented in 
Paper I: 



logiVcr = (2.82 ± 0.22) + (0.049 ± 0.004) x Fa 



(1) 



is obtained, with a correlation coefficient of 0.96 and a = 
0.46. Several data points in Fig. 5 (left) correspond to lines 
not included in the HR grating wavelength ranges discussed 
here, having been obtained from the lower resolution (LR) 
grating spectra of Paper I (including [S ll] AA6717, 6730, 
[CI III], [Ar III], [Ne iii], and [N ii] A5755). In all cases the 
measurements are for the rest velocity component of the 
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Figure 2. Spectra of the background subtracted LV2 emission (black) and local M42 vicinity (red) from various high resolution VLT 
FLAMES gratings smoothed with a three point average filter for illustration purposes (flux brightness per 0.31x0.31 arcsec^). Prominent 
lines have been identified in blue for lines common to both LV 2 and M42, in red for lines dominated by M42 emission. Green symbols 
mark lines arising from LV 2's bipolar jet. 
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Figure 3. As in previous figure. 
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Figure 4. A comparison of line profiles from the intrinsic LV2 spectrum (black) and from the local Orion nebula (red). The peak flux 
has been scaled to unity. On the LV2 spectrum note the extreme weakness of [Ar iv] A4711 and of the [O ll] AA3726, 3729 lines due to 
coUisional de-excitation at high densities, as well as the relative strength of the three [Fe m] A4658 components. 
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Figure 5. (Left): Plot between the critical density (at 8000 K) of various emission lines versus the fraction (in per cent) of specific line 
fiux emitted from LV 2 (background-subtracted spectrum) over that emitted from LV 2 + M42 (observed spectrum) . The red data points 
(crossed-circles) are from the HST FOS measurements of Paper I and have been excluded from the linear fit (dashed line); the rest are 
from the FLAMES/Argus observations. (Right): The excitation temperature of a given line's upper level is plotted against the same 
quantity and for the same lines as in the left panel. 



lines (excluding the outlying jet lobes which correspond to 
components 1 and 4 in Table 2). No correlation is evident 
when Fa is plotted against the lines' excitation temperature 
(Fig. 5 right) ; this is to be expected as LV 2 if of similar Tc 
as the local nebula (within ft;800K). It is therefore estab- 
lished that due to its high electron density (~10® cm~^), 
the surface brightness of the proplyd compared to the local 
Orion nebula is very high in lines with large critical densi- 
ties such as the auroral [O ill] and [N ll] lines, and the UV 
transitions. For the O ll A4649 and C ll A4267 lines F\ = 
71 and 68 per cent. This is essentially identical to the value 
found for H/3 and other Hi Balmer lines. The significance of 
this result is discussed in the conclusions. 



3.2 Line maps, velocities and widths 

The profiles of lines emitted from the proplyd as opposed 
to those emitted from the local Orion nebula are compared 
in Fig. 4 in the heliocentric velocity frame. Gaussian fits to 
the profiles have been employed to measure the line helio- 
centric velocity and FWHM for the intrinsic proplyd spectra 
and for the co-added background M42 spectrum (Table 2). 
Regarding the intrinsic LV 2 spectrum, double Gaussian pro- 
files were fitted to lines of high S/N ratio (components 2 and 
3 in Table 2 corresponding to the core or 'cusp'), while single 
Gaussians were fitted to the core of weaker lines. The out- 
lying blue- and red-shifted jet lobes were each fitted with 
a single Gaussian (components 1 and 4 respectively in Ta- 
ble 2). Lines such as [N ll] and [Ar iv] whose profiles are 
double-peaked lack FWHM measurements, and their helio- 



centric velocities correspond to that of their central troughs. 
The fit parameters pertaining to the core of [O ill] A4959 are 
less accurate due to the increased uncertainty of the exact 
line profile in the 20-40 kms~^ interval (marked by ':' in 
Table 2). The FWHMs in Table 2 have been corrected for 
instrumental broadening by quadrature subtraction, but the 
thermal (Doppler) broadening has not been corrected for: at 
ion temperatures of 9000 K the thermal width of the H I and 
He I lines is 20 and 10 kms~^ respectively. The hyperfine 
structure causes additional broadening to lines such as C ll 
A4267. 

The widths of lines from LV 2 are typically larger than 
those from the nebula. This is an expected result (e.g. Hen- 
ney & O'Dell 1999) and is due to the fact that the proplyd 
is experiencing continuous mass-loss through a photoevap- 
orated fiow arising from an embedded circumstellar disk by 
the action of far-UV radiation (e.g. Johnstone, HoUenbach & 
Bally 1998). The fiow is then converted into a photoionized 
wind by the Lyman continuum extreme UV (EUV) radia- 
tion field of 8^ OriC. The ionized wind and any externally 
irradiated ionized jets, when present as in LV2, emit essen- 
tially an H II region spectrum (see Figs. 2, 3). The [O ll], 
[N ll] , and [Ar iv] lines show the efi'ects of coUisional quench- 
ing as discussed above. It is noteworthy that the C ll and 
O II recombination lines arising from LV 2 appear very sim- 
ilar to each other and to the [O in] coUisional lines. In the 
case of O ll versus [O ill] this could mean that the lines are 
emitted from gas volumes of rather similar physical prop- 
erties. For example, this could be linked to the fact that 
both sets of lines return similar abundances for O^"*", that 
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Figure 6. Heliocentric velocity of various lines emitted from ions 
of lower to higher degree of ionization (left to right on the horizon- 
tal axis; the hydrogen lines are all due to recombining H+). The 
data for LV 2 (black colour) were obtained from the background- 
subtracted spectrum, while those for M42 (red colour) are from 
the co-added spectrum as defined in the text. For [N ll] and [O in] 
the symbols at higher velocity correspond to the auroral lines 
(A5755 and A4363 respectively), while those at lower velocity cor- 
respond to the nebular lines (A6584 and A4959). [N ll] A5755, [O l] 
A6300 and [S III] A6312 were measured on the LR5 grating spectra 
of Paper I. 



is, they yield a low abundance discrepancy factor (ADF ~ 1, 
defined as the ratio of the recombination Une (RL) over col- 
Usional line (CL)-derived O'^^ abundance; see Paper I). In 
contrast, there have been cases of planetary nebulae such as 
NGC 6153 and NGC 7009 where the smaller widths of O ll 
RLs versus [O ill] CLs have been linked to the presence of 
cool, H-poor plasma and high ADFs for O'^"'" (e.g. Barlow et 
al. 2006). 

The velocities and FWHMs of O ll RLs and [O in] CLs 
in the M42 spectrum are also very similar. The Doppler- 
corrected widths of H I and metal lines in M42 are compa- 
rable. There is some indication for slightly larger Doppler- 
corrected linewidths of Hi compared to [S ll] or [O ll], but 
not when compared to [O ill]. In an extensive study of the 
gas velocity structure in M42 at an effective resolution of 10 
kms~^, Garci'a-Diaz et al. (2008) concluded that the non- 
thermal broadening is not systematically different between 
recombination lines and coUisional lines. It is also notable 
that the velocities of the [O ll] violet doublet components 
are separated by 3-4 kms~^ in the M42 spectrum; the same 
is true for the [O ill] nebular and auroral transitions. This 
points to a degree of variation in the ion density and/or tem- 
perature along the line of sight which warrants further inves- 
tigation. In Fig. 6 the behaviour of the radial velocity of lines 
emitted from low to highly ionized species is examined for 
LV 2 and for the adjacent M42 nebular field. A shallow trend 
of decreasing velocity for increasing ionization degree is seen 
for LV 2, but the [S ii] A4069, [O iii] A4363 and [S iii] A6312 
data which are due to lines of high critical density muddle 
the picture. The velocities plotted are those of component 2 



(Table 2). The trend does not strengthen if instead we were 
to adopt velocities obtained from single component Gaus- 
sian fits to the cusp of LV2 (i.e., treating components 2 & 3 
as one component), but in that case LV 2 appears red-shifted 
versus M42 for all lines (by 6 ± 3 kms~^). The scatter of 
the data is larger for M42 (even amongst the H I transitions 
which are emitted from the H+ recombining ion), and this 
is probably real as the observational uncertainties are typ- 
ically much smaller. The behaviour of ion velocity versus 
ionization potential (IP) for proplyds and their environs has 
been explored by Henney & O'Dell (1999) who found trends 
of increasing blue-shift for increasing IP in three proplyds, 
while the opposite trend was observed in a fourth proplyd. 
Regarding M42, the trend of decreasing velocity for higher 
IP is commonly explained in terms of the ionization stratifi- 
cation of the nebula with low ionization lines being emitted 
closer to the ionization front (IF) which lies in front of the 
molecular cloud, and highly ionized species existing closer 
to the ionizing star having been accelerated away from the 
IF towards the observer (e.g. Esteban & Peimbert 1999 and 
references therein). 

In Fig. 7 line maps of LV 2 and its vicinity are shown, 
created from the Ha, H7, and [Fe ill] A4658 surface bright- 
ness distributions in selected velocity ranges. For the H I 
maps, sections of the field clear of proplyd emission were 
fitted with Legendre polynomials and the resulting model 
surfaces were subtracted from the original frames to yield 
background-subtracted images. The emission from the red- 
and blue-shifted lobes of LV 2's jet is marginally spatially 
resolved: a separation of ~1.5 spaxels (0.47 arcsec) is mea- 
sured between the respective emission peaks of the green and 
blue contours. The Ha channel maps were obtained from the 
highest velocity resolution spectra of this data set and show 
in colour the Vhoi = 0-60 and 60 - 110 kms~^ range in 
separate panels with the bipolar jet emission displayed as 
contours. The proplyd appears more diffuse and elongated 
towards the tail in the latter velocity interval which encom- 
passes less than 25 per cent of the line luminosity, as com- 
pared to a more compact shape close to the line centre. This 
would suggest that the ~60 - 110 kms^^ interval intermedi- 
ate between the cusp of LV 2 and the fast red-shifted jet lobe 
is occupied by gas which flows away from the cusp and goes 
on to form the proplyd's comet-like tail. This corroborates 
the study of Henney et al. (2002) who hinted at this pos- 
sibility. The same is also evident in the H7 map where the 
colourmap displays the whole - 110 kms~^ interval. From 
the Hi maps and a comparison between the shape of the jet 
contours with the orientation of the tail it would seem that 
the projected axis of the bipolar jet is at an angle of «90 
deg with the tail of LV 2. 



3.3 Physical conditions and jet properties 

The red-shifted lobe of LV2's jet is prominent in several 
emission lines from hydrogen, helium and heavier species 
(0°, N+, S+, S^+, 0^+, Ne^+, Fe^+). The blue-shifted lobe 
is also well-detected in [O i], [O iii], [S iii], [Fe iii]A4658, and 
the strong Hi lines (see Figs. 2, 3 of this paper and fig. 7 
of Paper I): the lines are however a factor of a few weaker. 
A compilation of line fiuxes from the high-dispersion HR 
spectra of this paper and the lower-dispersion LR spectra of 
Paper I is presented in Table 3 for the red-shifted lobe. 
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[Fe HI] 4658 A (16, 45) km/s 




xlO 
2.5 
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Figure 7. Maps of LV2 and vicinity in Hi and doubly ionized iron emission. Reddening correction has not been applied. The colourmap 
units are ergs~^ crn"^ per 0.31" X 0.31" spaxel for the velocity range noted at the top of each panel. The various contours represent 
emission from the proplyd's bipolar jet in selected heliocentric velocity ranges (see bottom of each panel). The contours in both Ho panels 
are at the 5.0x10-^5^ 2.0x10-" (green) and 2.0x10-^5, 5.0x10-^5 (blue) levels. The H7 contours are at the 2.0x10-1-*, LOxlQ-^^ 
(green), and 1.5x10-", 3.0x10-" (blue) levels. The [Fe III] contours are at the LOxlO-l^ (blue) and LOxlQ-", 3.0x10-" (green) 
levels. The H I maps have been background-subtracted. For display purposes the maps have been reduced from their original size of 
14x22 spaxcls. 
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Table 2. Line structure in the intrinsic LV2 spectrum as compared to the local M42 field." 











LV 2 proplyd 








Component 


1 


2 


3 


4 


1 


2 


3 


4 






Vq (kms-i) 






FWHM 


(kms-i) 




H12 3750 


- 


28.7±0.4 


- 


- 


- 


41.2±0.9 


- 


- 


Hll 3770 


- 


28.4±0.3 


- 


- 


- 


44.6±0.8 


- 


- 


HIO 3798 


- 


28.9±0.3 


- 


- 


- 


47.2±1.0 


- 


- 


H9 3835 


- 


28.1±0.5 


- 


- 


- 


46.1±0.7 


- 


- 


[S II] 4069 


- 


28.8±0.3 


- 


129.8±1.4 


- 


24.0±0.8 


- 


52.7±4.8 


[S II] 4076 


- 


28.7±0.5 


- 


- 


- 


36.2±0.8 


- 


- 


US 


- 


19.0±0.5 


46.1±0.6 


133.9±0.5 


- 


34.4±0.8 


34.8±0.8 


48.7±1.2 


C II 4267 


- 


34.8±0.9 


- 


- 


- 


52.6±2.0 


- 


- 


H7 


-105.6±0.6 


23.6±0.3 


41.6±1.3 


141.8±0.5 


62.1±1.7 


30.9±1.2 


45.1±0.9 


47.4±1.0 


[O III] 4363 


-112.7±5.5 


24.7±0.5 


41.9±2.8 


143.8±0.5 


73.9±14.4 


22.3±1.8 


36.3±2.6 


38.6±1.1 


[Fe III] 4607 


- 


- 


- 


150.8±1.0 


- 


- 


- 


37.1±2.5 


O II 4649 


- 


20.6±0.6 


54.2±1.1 


- 


- 


28.0±2.7 


34.5±5.4 


- 


O II 4650 


- 


32.3±3.2 


- 


- 


- 


69.0±8.0 


- 


- 


[Fe III] 4658 


-106.9±2.1 


26.4±1.7 


- 


144.9±0.2 


58.6±5.0 


25.6±3.9 


- 


33.1±0.5 


O II 4661 


- 


32.8±2.4 


- 


- 


- 


62.2±6.0 


- 


- 


[Fe III] 4702 


-114.6±2.0 


17.5±2.0 


- 


139.4±0.2 


69.0±5.3 


25.7±4.7 


- 


32.0±0.6 


He I 4713 


- 


21.3±0.6 


39.7±3.9 


- 


- 


24.9±2.9 


41.6±3.5 


- 


[Ar iv]4740 


- 


35.2: 


- 


- 


- 


- 


- 


- 


He I 4922 


- 


30.3±0.5 


- 


- 


- 


46.9±1.3 


- 


- 


[0 III] 4959 


-109.0±2.0 


25: 


55: 


138.8±0.5 


45.2±2.3 


36: 


28: 


42.1±1.1 


[Fe III] 5011 


- 


- 


- 


142.1±0.5 


- 


- 


- 


35.8±1.1 


[N II] 6584 


- 


39: 


- 


139.0±1.0 


- 


- 


- 


39.6±2.0 


Ha 


-102.9±2.2 


27.0±2.1 


43.2±7.4 


142.7±0.5 


51.9±1.5 


36.7±6.7 


47.5±3.6 


50.0±3.0 










Local M42 nebula 












Vq (kms-l) 






FWHM 


(kms-l) 




[O II] 3727 


- 


25.0±0.2 


- 


- 


- 


14.3±0.3 


- 


- 


[O II] 3729 


- 


21.8±0.3 


- 


- 


- 


15.7±0.4 


- 


- 


H12 3750 


- 


21.0±0.3 


- 


- 


- 


25.8±0.3 


- 


- 


Hll 3770 


- 


21.0±0.2 


- 


- 


- 


26.0±0.5 


- 


- 


HIO 3798 


- 


20.3±0.4 


- 


- 


- 


26.5±0.5 


- 


- 


H9 3835 


~ 


20.3±0.3 


- 


- 


- 


25.4±0.4 


- 


- 


[S II] 4069 


- 


25.4±0.3 


- 


- 


- 


12.5±0.6 


- 


- 


[S II] 4076 


- 


25.4±0.5 


- 


- 


- 


11.5±0.8 


- 


- 


m 


- 


21.8±0.1 


- 


- 


- 


25.4±0.4 


- 


- 


C II 4267 


- 


29.3±0.5 


- 


- 


- 


27.8±1.3 


- 


- 


H7 


- 


27.7±0.0 


- 


- 


- 


26.0±0.1 


- 


- 


[O III] 4363 


- 


28.4±0.2 


- 


- 


- 


15.3±0.5 


- 


- 


[Fe III] 4607 


- 


34.6±0.6 


- 


- 


- 


16.9±1.4 


- 


- 


[O II] 4649 


- 


26.8±0.2 


- 


- 


- 


16.6±0.5 


- 


- 


[O II] 4650 


- 


26.6±0.4 


- 


- 


- 


18.6±1.0 


- 


- 


[Fe III] 4658 


- 


29.1±0.2 


- 


- 


- 


15.9±0.5 


- 


- 


O II 4661 


- 


24.9±0.6 


- 


- 


- 


18.6±1.4 


- 


- 


[Fe III] 4702 


- 


25.5±0.2 


- 


- 


- 


14.6±0.6 


- 


- 


He I 4713 


- 


25.1±0.2 


- 


- 


- 


20.1±0.3 


- 


- 


[Ar iv] 4740 


- 


30.4±0.3 


- 


- 


- 


14.9±0.7 


- 


- 


He I 4922 


- 


22.7±0.1 


- 


- 


- 


18.9±0.3 


- 


- 


[O III] 4959 


- 


24.3±0.2 


- 


- 


- 


16.2±0.4 


- 


- 


[Fe III] 5011 


- 


26.1±0.4 


- 


- 


- 


12.0±0.9 


- 


- 


[N II] 6584 


- 


28.3±0.2 


- 


- 


- 


11.1±0.4 


- 


- 


Ha 


- 


29.7±0.1 


- 


- 


- 


26.5±0.3 


- 


- 



" Velocities are heliocentric. The FWHMs have been corrected for instrumental broadening using the Th-Ar arc line measurements 
reported in Section 2. Highly uncertain values arc marked with ':'. 



In order to compensate for the fact that there are gaps 
in the wavelength coverage of the FLAMES HR gratings 
and that the LR and HR observations were taken on differ- 
ent nights and under different observing conditions, the HR 
measurements have been cast in units relative to H/3 using 
the LR spectra of Paper I in the following manner: 



f(A)jct r f (A)jot 

F(H/3)jet Li^(A)corc 



r f (A)corc 1 
Hr/ LF(H/?)jetJLR; 



(2) 



The [Fe ill] A4658 line was scaled via He I A4713, and [Fe ill] 
A5011 was scaled via He I A5016. Ha was scaled via [N ll] 
A6584 as the rest velocity component of Hq probably suf- 
fers from saturation effects on the HR14 and LR5/6 grating 
spectra. The relative line fluxes were then dereddened with 
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c(H^) = 1.2±0.2 using the modified Cardelli, Clayton & 
Mathis (1989; CCM) law from Blagrave et al. (2007), with a 
total to selective extinction ratio Rv = 5.5. This reddening 
constant was obtained from the H(5/H/3 and H7/H;3 ratios 
from a comparison of the observed ratios to their theoretical 
values from Storey & Hummer (1995) at Tc = lO"* K and N^ 
= 10® cm"''. The observed Ha/H/3 ratio is considered unre- 
liable due to Ha likely being saturated. Using the same H I 
line ratios, c(H/3) for the blue-shifted jet lobe is measured to 
be 0.9 and 1.5 respectively, in approximate agreement with 
the reddening deduced for the red-shifted lobe above and for 
the cusp of LV 2 from Paper I. 

In Fig. 8, a group of curves is shown on the (To, A^o) 
diagnostic plane representing solutions for the dereddened 
line ratios from Table 3 applicable to the red-shifted jet. 
The ratios [N ii] A5755/A6584, [O i] A5577/A6300, [O iii] 
A4363/A5007, and [Fe iii] A5011/A4658 (^Pi-^D2)/(^F4- 
^D4) were used, employing the multi-level ion modelling pro- 
gram EQUIB developed at University College London. The 
C III] A1907/A1909 curve is from the HST STIS observa- 
tions of LV2 by Henney et al. (2002). From Fig 8 we de- 
duce that the mean electron temperature and density of the 
jet are 9000 - 10,000 K and ~10'' cm"^, respectively. The 
iron lines indicate that overall higher densities of 2.5x10® 
cm""' are appropriate for their emitting region, and densi- 



ties log (A'c/cm 



7.0_o g are obtained from the [Fe iii] 



A4702/A4658 (■'Fg- D3)/(^F4- D4) ratio when comparing 
with the theoretical ratios of Keenan et al. (2001). The ion- 
ized carbon, nitrogen, and neutral oxygen lines also indicate 
densities of ~10® cm""'. The fact that a range of densities 
are obtained from various diagnostics is an indication that 
the jet is stratified in density and probably has a clumpy 
small-scale structure. This is corroborated by the presence 
of neutral oxygen in the body of the jet, already noted in Pa- 
per I (see fig. 7), indicating that dense quasi-neutral clumps 
can withstand complete photoionization in the jet fiow. 

Our analysis suggests that previous mass-loss rate esti- 
mates from the jet of ~10~* M0 yr~^, which were based on 
considerations of a fully ionized flow (e.g. Henney et al. 2002; 
Vasconcelos et al. 2005), would actually be lower limits. 



3.3.1 Velocity-resolved physical conditions 

We used the high dispersion [O in] line profiles obtained with 
the VLT Argus fibre array in conjunction with the C ill] HST 
STIS observations by Henney et al. (2002) from programme 
8120 to map the variation of electron temperature and den- 
sity as a function of the heliocentric velocity of the ionized 
gas in LV2. These C ill] and [O ill] lines should originate 
from roughly the same volume of highly ionized gas. The re- 
sults from this analysis are shown in Fig. 9 where the [O in] 
profiles are from the intrinsic (background-subtracted) LV 2 
spectrum. The [O ill] and C ill] profiles were first binned 
to a velocity resolution of 3 kms~^ pix~^ and the electron 
density was computed as a function of Viiei from the density- 
sensitive (reddening and Te-independent) C ill] A1907/A1909 
ratio for an electron temperature of 9000 K. The obtained 
density profile (Fig. 9 middle panel) was then used to com- 
pute the temperature profile from the [O ill] auroral to neb- 
ular line ratio (Fig. 9 top panel). The density variations that 
are observed across the bulk of LV 2's emission in the 10-80 
kms~^ range are essentially the same as those obtained by 



Henney et al. (2002) who analyzed the C ill] data in a sim- 
ilar way. The electron density in LV2 peaks very close to 
the rest velocity of the [O ill] lines [logA'c (cm"') = 6], then 
falls off towards a local minimum at ~100 kms"^, before 
rising again to values of log [A^e (cm^)] = 6 and higher in 
the velocity range occupied by the red-shifted jet emission. 

The analysis reveals that steep temperature variations 
are observed: (i) in the velocity interval between the bulk 
LV 2 emission and the onset of the red-shifted jet lobe, (ii) 
at higher velocities beyond the 'red' edge of the jet. The ob- 
served amplitude is ~3500K with a maximum of 12 500 K 
on the 'blue' side of the jet, while even higher tempera- 
tures (admittedly with very large associated uncertainties) 
are observed on its 'red' side. These steep Tc variations could 
be the signature of a shock discontinuity and its associated 
heating/cooling zone which is kinematically resolved here. 
A shock discontinuity can form at the locus where the jet 
encounters and bursts out from the main ionization front 
surrounding LV2 (the H'^/H"'' contact surface) and/or at 
the outer physical envelope of the conical/cylindrical jet. 

On the other hand, at velocities associated with the 
bulk emission near the rest frame of LV 2, and with the bulk 
jet emission, the temperature is ~9000 K which is typical of 
photoionized gas. This measurement (the horizontal black 
line in Fig. 9 top) is corroborated by the analysis of the in- 
tegrated line ratios for the rest velocity component (9000 ± 
600 K) from the lower-dispersion spectra of Paper I (table 
3), and from the diagnostic diagram pertaining to the jet in 
the present work (Fig. 8). In the temperature plot of Fig. 9 
the red data points correspond to the intrinsic [O ill] line 
profiles where the interval 20-40 kms~^ has been excluded 
from the analysis due to the increased uncertainty in the ex- 
act profile of A4959 very close to the bulk velocity of M42. 
Specifically, Tc = 9240 with cr = 815K for V^a £ [0, 80] 
kms~^ across LV2's cusp (where the scatter indicates real 
variations rather than observational uncertainties). This ve- 
locity interval contains >98 per cent of the intensity of the 
[O III] auroral and nebular lines near the rest frame of LV2, 
excluding the outlying jet lobes. The analysis was also per- 
formed using the observed (LV 2 -I- M42) line profiles re- 
sulting in a variation (black crosses) that follows closely the 
previous one across most of the velocity range sampled. In 
this case we find that across the cusp Tc = 9020 with a = 
910 K for VhoI £ [0, 80] kms"\ We therefore deduce that 
the bulk of the [O ill] emitting region in LV2 appears to 
be fairly isothermal with a <10 per cent rms temperature 
variation, while considering the uncertainties the main body 
of the jet may be less so. 

Using the information obtained above for the physi- 
cal conditions in the red-shifted jet lobe we have computed 
ionic abundances for gas temperatures of 9000 and 10 000 K, 
adopting an electron density of 10® cm""'. The results are 
presented in Table 4. Both dominant ionization states of sul- 
phur are observed allowing us to obtain its gas-phase abun- 
dance. We have assumed an ionization correction factor {icf) 
of 1.05 meaning that higher ionization stages probably con- 
tribute less than five per cent of the total (eq. 5, Paper I). 
This icf is applicable to the cusp of LV2 and can be consid- 
ered to be a valid assumption for the jet also given the very 
similar ion temperature and density in the two. The abun- 
dance of sulphur is then 12-|-log(S/H) = 6.93 ± 0.06, that 
is, ~0.2 dex less than in the Sun. This is in good agreement 
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Figure 8. Electron temperature and density solutions for the 
red-shifted lobe of LV2's jet. Dotted lines bracketing thick solid 
lines of the same colour correspond to Icr min/niax values of the 
respective diagnostic ratio. The C III] curve is from the HST STIS 
data of Henney et al. (2002) discussed in the text. 



with the abundance of 6.83 measured from the rest velocity 
component of the lines due to emission from the cusp (table 
7, Paper I). An independent determination for the gas-phase 
S/H ratio in M42 is 6.89±0.02 based on SST data (Rubin 
et al. 2011). 



3.3.2 The Fe abundance and dust considerations 

Of particular interest is the abundance of Fe'^"''/H^ in the jet 
which places a lower limit to the total iron abundance. We 
used a 25-level ionic model for Fe ill employing the transition 
probabilities of Nahar & Pradhan (1996) and the collision 
strengths of Zhang (1996) to compute the abundance of Fe^^ 
using the AA4658, 5011 lines, finding 12 + log(Fe^+/H+) = 
7.01±0.05. This is 2 dex higher than the total Fe/H abun- 
dance in the cusp of LV 2 (Paper I). It is therefore established 
that (i) iron in the main body of LV 2 is depleted by a factor 
of at least 100 with respect to that in the jet, and (ii) there 
is a factor of at least nine more iron in the jet (>1 dex) than 
in the bulk M42 gas (cf. Paper I and Esteban et al. 2004). 
The jet therefore has an iron abundance of >30 per cent 
solar or a depletion relative to solar [Fe/H] = log(Fe/H) — 
log(Fe/H)Q > —0.5 dex. If we adopt the same icf(Fe) for 
the jet as that for the core from Paper I then [Fe/H]jot '^ 
-0.2 dex. 

An increase in the gaseous Fe abundance associated 
with the high-velocity HH 202 outflow in the Orion neb- 
ula has also been attributed to dust destruction by Mesa- 
Delgado et al. (2009), who found —(0.3-0.5) dex of iron de- 
pletion with respect to solar. The depletion of iron in the 
core of LV2 sampled by the rest velocity line component 
is [Fe/H]coro = -(2.54±0.40) dex (Paper I). This is higher 
than the depletions of (i) —1.3 dex associated with the low- 
velocity component of the HH 202 flow (Mesa-Delgado et al. 
2009), and (ii) —1.5 dex measured in the local background 
of LV 2 which samples the generic M42 gas (Paper I) . The 
iron depletion in LV2's rest frame (0.003-0.007 of the solar 



..•-• 



03 




200 



Figure 9. The electron temperature (top) and electron density 
(middle) distribution in velocity space along the line of sight to- 
wards LV2, based on background-subtracted [O III] and C III] 
observations (as shown in the bottom panel and discussed in the 
text respectively). The solid/dashed horizontal lines in the top 
panel denote the temperature and Icr uncertainty measured from 
the integrated [O III] line ratio for the cusp in Paper I (table 3). 



abundance) is also higher than in typical Galactic H ll re- 
gions, and lies within the high end of the range exhibited 
by some planetary nebulae (Delgado Inglada et al. 2009). In 
terms of the general interstellar medium this level of deple- 
tion is comparable only to what is found towards a minority 
of sight lines, including the direction to the star ( Oph which 
represents the prototype of strong depletions observed in the 
cold ISM (Jenkins 2009). 

The fact that sulphur, which is a volatile element [its 50 
per cent condensation temperature, Tc, into troilite (FeS) is 
664 K; Lodders 2003], shows essentially the same abundance 
in the core and jet of LV2, whereas the refractory species 
iron (Tc = 1334 K) is much more abundant in the jet than in 
the core leads to the following conclusions: (i) iron-bearing 
condensates must be very efficiently destroyed in the fast 
jet, probably as a result of grain sputtering processes, thus 
substantially enhancing the overall gaseous abundance of Fe 
(factor of X 220 increase over LV 2's rest frame, and x 19 over 
M42); (ii) there is no evidence for a significant depletion of 
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sulphur onto dust and therefore LV 2 has a genuinely lower 
S abundance than the Sun (—0.2 dex). 



Table 3. Fluxes of emission lines from the red-shifted jet of LV 2 
in units such that H/3 = 100." 



3.3.3 The Ne/S abundance ratio 

The Ne/S abundance ratio in LV2 is well constrained via 
the proxy ionic ratio Ne'^'''/S^''' as it involves lines of similar 
excitation energy (37.2 kK for [Ne ill] A3869 and 39.1 kK 
for [S in] A6312) and critical density (7.1x10*' vs. 1.0x10^ 
cm^'^), and is therefore fairly insensitive to the assumption 
of Tc or A'e for their emitting zones. This ratio is 13.5 and 
14.8 in the jet and core of LV 2 respectively (cf. Table 4 with 
Paper I, table 5) whereas the canonical solar Ne/S ratio is 
6.5 (Asplund et al. 2009). Our measurement for LV2 agrees 
with the value of 13.0 applicable to M42 from infrared lines 
measured by Spitzer (Rubin et al. 2011). The IR lines set 
a good benchmark as their are immune to Tc or reddening 
biases. It is discomfiting however that the ratio in the local 
vicinity of LV2 measured using the optical lines is only 4.9 
(Paper I, table 5), and in agreement with the value of 4.8 
from Esteban et al. (2004) applicable to a different M42 zone 
obtained also from the optical lines. 

Two factors could conspire to bring about a discrepancy 
between the optical and IR values of the Ne/S ratio for M42: 
in terms of Te effects, the offset could be explained if the op- 
tical Ne'^^/S'^^ ratio was underestimated due to a variation 
of ~1900K between a cooler Ne'^'*' and a warmer S^"*" emit- 
ting volume (rather than assuming the same temperature 
for both). Secondly, the differential reddening correction be- 
tween 3869 and 6312 A is a factor of about two, and if the 
extinction along the optical observation sightlines has been 
underestimated in any way this would also cause the optical 
Ne^"'"/S'^^ ratio to appear too low. On the other hand, it 
would be very hard to appeal to Tc (or extinction) biases 
to explain the opposite, i.e. that the Ne/S ratio in LV2 has 
been overestimated, as this would have to act opposite to 
the natural tendency of the A3869 and A6312 lines to be 
more efficiently excited at lower and higher To, respectively. 

We therefore conclude that the Ne/S ratio in LV2 is 
genuinely higher than solar by factors of 2.1 - 4.3. The lower 
value is obtained from the Ne^^/S^^ determination applica- 
ble to the jet and core of LV2. The upper limit is obtained 
from the total Ne/S ratio for the rest velocity component 
(core) from Paper I and may be affected by ionization cor- 
rection factor considerations. This adds to the body of ev- 
idence that the canonical solar Ne (Asplund et al. 2009) 
appears to be low compared to M42 (Rubin et al. 2011) and 
the LV2 proplyd (Paper I and this work) based on consid- 
erations of the Ne/S ratio, the absolute neon abundance in 
B-type stars in Orion (Simon-Diaz & Stasihska 2011 and 
references therein), or when comparing the solar Ne/O ra- 
tio to the mean value from Galactic planetary nebulae and 
H II regions (e.g. Wang & Liu 2008). A new solar abundance 
model by Delahaye et al. (2011) partially rectifies this situ- 
ation by raising the Ne abundance to 8.15±0.17 (-1-0.22 dex 
with respect to the value in Asplund et al 2009). 



4 DISCUSSION AND CONCLUSIONS 

We have presented spatially resolved, high dispersion spec- 
tra of LV2 and its Orion nebula vicinity based on VLT 



Line 



F{\) 



/(A) 



[Ne oi] A3869 
[S II] A4069 
H<5 A4101 
H7 A4340 
[O III] A4363 
He I A4471 
[Fe III] A4658 
[Fe III] A4702 
[Fe III] A4734 
H;9 A4861 
[O III] A4959 
[O III] A5007 
[Fe III] A5011 
[O i] A5577 
[N 11] A5755 
[O i] A6300 
[S III] A6312 
H« A6563 
[N 11] A6584 



45.1 ± 6.45 

5.15 ± 0.71 

20.4± 0.9 

32.8 ± 2.6 

5.93 ± 0.55 

5.21 ± 0.55 

12.4 ± 0.94 

5.56 ± 0.57 

3.03 ± 0.39 

100 ± 4 

103 ± 9.4 

309 ± 29 

4.75 ± 0.72 

0.231± 0.046 

6.97 ± 1.16 

9.10 ± 0.49 

6.99 ± 0.50 

287 ± 20: 

7.21 ± 0.43 



67.2 ± 9.7 

7.19 ± 0.99 
27.9± 1.6 
42. 3± 3.4 

7.44 ± 0.70 
6.22 ± 0.66 
13.0 ± 1.0 
5.9 ± 0.61 

3.20 ± 0.41 
100± 5 

99. 0± 9.0 

290± 28 

4.46± 0.68 

0.174± 0.035 

4.95± 0.83 

5.39± 0.28 

4.14± 0.30 

201± 14: 

3.90± 0.25 



" Several lines were measured on the lower resolution Argus 
spectra using the gratings quoted in table 1 of Paper I: [Ne III] 
A3869 (LRl), He i A4471 (LR2), H/3 A4861 (LR3), [O i] A5577 
(LR4), [N 11] A5755 (LR4/LR5), [O i] A6300 (LR5), [S iii] A6312 
(LR5). The values for Ha are unreliable. See the text for details. 



Table 4. Abundances for the red-shifted jet lobe of LV2 (in a 
scale where logH = 12 for the first seven rows and in decimal 
units for the last three rows)." 



Species 



To = 9000 K 
A^o = lO** cm-3 



To = 10000 K 
Ac = 10^ cm-3 



N+/H+ 

02+ /H+ 

Ne2+/H+ 

Fe2+/H+ 

S+/H+ 

S2+/H+ 

S/H 

Fe2+/02+ 

Ne2+/02+ 

Ne2+/S2+ 



8.60"' 



7 7Q-I-0.07 
'•"J-0.08 
^+0.04 
-0.05 
o nc;+0.10 
«-U3-0.03 
a+0.03 
-0.03 
.+0.05 



7.08"* 



6.24 
6.92 
7.02 



0.07 
+0.04 

0.04 
+0.04 

0.04 



7 46+0-°'' 
'•*°-0.08 
3+0.03 

-0.05 
-+0.10 

-0.03 
,+0.04 

-0.03 
7+0.08 

-0.06 
1+0.05 

-0.03 
,+0.05 



8.29: 
7.85^ 
6.92^ 
6.071 
6.71^ 



6.82] 



0.05 



0.030±0.004 
0.28±0.08 
13.5±3.7 



0.043±0.006 
0.36±0.10 
13.8±3.8 



" Solar abundances for comparison: N = 7.83±0.05, O = 
8.69±0.05, Ne = 7.93±0.10, S = 7.12±0.03, and Fe = 7.50±0.04 
(Asplund et al. 2009). 

observations. The deep spectra enabled the detection and 
resolution of emission lines arising from the proplyd's irra- 
diated bipolar microjet. Maps of LV 2 in the light of H l and 
[Fe III] lines have been presented, where the jet lobes appear 
spatially resolved from the core of the proplyd. 

An abundance analysis was undertaken for the brighter 
red-shifted jet lobe complementing the abundance analysis 
of the core of LV2 presented in Paper I. A noteworthy re- 
sult is the significantly enhanced abundance of iron in the 
jet versus the core of LV 2, pointing to very efficient dust de- 
struction mechanisms in fast irradiated microjets. The over- 
all depletion of iron in the proplyd's rest frame is higher 
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than in Galactic H ll regions and is comparable to what is 
found along cold ISM sight lines. The Ne/S abundance ra- 
tio in LV2 as probed by optical lines is >2 times solar, in 
agreement with the value applicable to the generic M42 gas 
from IR line measurements by Spitzer. 

The extreme-UV irradiated jet shows a range of ion- 
ization conditions allowing the survival of neutral oxygen 
pockets mixed in with highly ionized gas. A velocity-resolved 
analysis was performed probing the electron temperature 
and density across LV 2's rest frame and red-shifted jet lobe. 
The [O III] and C ill] emitting regions in LV 2 show a range of 
densities with a peak of l(f cm~^ near the line systemic ve- 
locity and an upper limit of ~10^ cm~^ from [Fe ill] emission 
within the red-shifted jet lobe. The electron temperature in 
the highly ionized zone in the main body of the proplyd 
shows a small 10 per cent rms variation. Larger tempera- 
ture variations associated with the onset of the jet emission 
have been attributed to a shock discontinuity. 

The far-UV to optical coUisional lines emitted from LV 2 
are correlatively quenched according to the critical densities 
of their upper levels, with lines of low critical density be- 
ing affected the most. This reinforces the conclusion of Pa- 
per I that in H ll regions containing dense zones, such as 
proplyds, filaments, globules etc., line ratios such as [O ill] 
A4363/(A4959-HA5007) and [N ii] A5755/(A6548-hA6584) wiU 
be heavily weighted towards the dense component of the gas. 
These line ratios are the main temperature diagnostics for 
Galactic and extragalactic H ll regions, whose validity rests 
on an independent way of determining the mean density of 
the plasma. This is usually done via the [O ll] or [S ll] dou- 
blet ratios which are demonstrably biased towards low den- 
sity regions (Fig. 5 left). The [N ll] ratio will be especially 
biased as the critical densities of its constituent lines are at 
a ratio of approximately 200 versus only 40 for the [O ill] 
lines, resulting in more uneven fractions of [N ll] A5755 ver- 
sus [N ll] A6584 to be emitted from dense and diffuse gas 
respectively, than is the case for [O ill] A4363 and [O ill] 
A5007 - this is demonstrated in Fig. 5. 

In this situation, and especially in cases of distant neb- 
ulae whose gas density distribution appears deceptively ho- 
mogeneous on large scales, or in the presence of conden- 
sations just below the angular resolution limit in nearby 
nebulae, the end result would be an overestimation of the 
electron temperature and an inevitable underestimation of 
the gas metallicity (e.g. Rubin 1989). Given the behaviour 
of the [N ll] and [O in] lines seen in Fig. 5 the former tem- 
perature would be more adversely affected that the latter in 
H II regions with a clumpy small-scale density distribution. 
It may be no coincidence therefore that Te([N ll]) values 
in H II regions are typically higher than re([0 ill]) values 
by ~1000K on average (e.g. Rodriguez & Delgado-Inglada 
2011). This is usually justified in the literature as being due 
to the fact that the low ionization zones in a nebula tend to 
have a higher temperature than the O ill zone. Based, how- 
ever, on the arguments presented here this behaviour could 
be at least partly spurious if the underlying density distri- 
bution of a nebula is sufficiently clumpy at sub-arcsecond 
scales. 

Regarding the abundance anomaly pertaining to RL 
versus CL diagnostics (e.g. Tsamis et al. 2003; Tsamis & 
Pequignot 2005; Garcia-Rojas & Esteban 2007; Ercolano 
2009; Esteban et al. 2009; Mesa-Delgado et al. 2010; Simon- 



Diaz & Stasihska 2011; Paper I), it should be noted that 
the fraction of O ll A4649, C ll A4267 and Hi RLs that are 
emitted from LV 2 is the same for all these lines (Ri70 per 
cent), while the rest is emitted from diffuse M42 gas along 
the line of sight. As a result the proplyd does not show el- 
evated ©^"""/H^ and C^^/H"'' abundance ratios, measured 
using RLs, compared to the Orion background (cf. Paper I). 
In our view this analysis goes some way towards explaining 
the abundance anomaly in H ll regions as a problem affected 
by the use and interpretation of the [O ill] and [N ll] plasma 
thermometers when applied to average spectra of large neb- 
ular volumes. 
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